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Abstract 

The pair-specific ground state energy eg{N) := E,g{N)/(N(N — 1)) of 
Newtonian N body systems grows monotonically in N. This furnishes 
a whole family of simple new tests for minimality of putative ground 
state energies E,g{N) obtained through computer experiments. Inspec- 
tion of several publicly available lists of such computer-experimentally 
obtained putative ground state energies Eg (N) has yielded several dozen 
instances of £-g{N) which failed one of these tests; i.e., for those N one 
concludes that £p(A) > £g(A) strictly. Although the correct E.g{N) 
is not revealed by this method, it does yield a better upper bound on 
&g{N) than 8,g{N) whenever £.g{N) fails a monotonicity test. The sur- 
veyed A^-body systems include in particular A'^ point charges with 2- 
or 3-dimensional Coulomb pair interactions, placed either on the unit 
2-sphere or on a 2-torus (a.k.a. Thomson, Fekete, or Riesz problems). 
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1 Introduction 

The pair-specific ground state energy eg{N) of Newtonian body systems 
exhibits the following monotonic dependence on A^: 

Proposition 1. Let A G M.'^ be a bounded and connected domain, and let 
Qj. G A. Assume the following hypotheses on U\{q, q): 

(HI) Symmetry: UAiq, q) = UA{q, q) 

{H2) Lower Semi- Continuity: UA{q,q) is l.s.c. on Ax A. 

For N > 2, define the pair-specific ground state energy by 

Then the sequence N £g{N) so defined is monotonic increasing. 

For the convenience of the reader we here reproduce the elementary proof 
from Appendix A in |Kie09j . 

Proof of Proposition [I].- 

We begin by noting that under hypotheses (-f^l) and {H2) the pair-specific 
ground state energy eg{N) defined in ([T]) is well-defined; i.e. eg{N) e M. 

To prove the monotonicity of A^ ^— > eg{N), with N > 2, let 8,g{N) denote 
the A^-body ground state energy, i.e. eg{N) = E,g{N)/ {N{N — 1)). Using the 
definition of 8,g{N) and the elementary graph-theoretical identity that the sum 
over all bonds in a complete A^-graph (A^ > 2) equals (A^ — 2)~^ the sum over 
all bonds of all its complete A^ — 1-subgraphs, and using the single inequahty 
that the minimum of a sum is not less than the sum of the minima, we find 



&g{N + l) = min V E 
^''^•■••'''~+i^<i<i<iv+i 



l<i<j<N+l 



l<i<j<JV+l 



>7^ E 

l<fc<iV+l 



{qi,...,<J^+J\{qJ^.|^^ 



N+l 

N- 



\EgiN). (2) 



Dividing by (A^ + 1)A^ yields ^^(A^ + 1) > eg{N), and the proof of the 
monotonicity of A^ ^— eg{N) is complete. □ 
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As already remarked in Appendix A of |Kie09j . Proposition [T] and its proof 
are quite elementary and presumably known, yet after a serious search in 
the pertinent literature I came up empty-handed, and additional consultation 
with several of my local expert colleagues have given me the impression that 
Proposition [1] is perhaps not known, and in any event not widely known. 

In this brief note we will be concerned with a very practical application 
of Proposition [1] which certainly is not generally known, as we will demon- 
strate. Namely, the monotonic increase with of the true pair-specific ground 
state energies Bg{N) furnishes a whole family of necessary criteria for mini- 
mality which any empirical list of computer-experimental data {£^(A^);A^ = 
2, 3, 4, A'^,} for such ground state energies needs to satisfy; put differently, 
we have the following sufficient criterion for failure to be minimal: 

VA^>2:(3n>l: e^g{N + n) < e^'giN) ^ TgiN) > 8.g{N)). (3) 

For each n > 1 one can use ([3]) as a test for any computer-experimentally 
produced list of putative ground state energies {£^(A^);A^ = 2, 3, 4, A'^,}. 
Of course, if any particular S,g{N) passes the n-th test for each n < N^: — N, 
i.e. if £g{N + n) > eg{N)\fn = 1, A",, — A^, this does not mean that this 
&g{N) is a true ground state energy; only failing a test is significant, for ([3]) 
asserts that a test-failing ^g{N) is surely not a ground state energy. Moreover, 
whenever ^g{N) fails any such test, then by Proposition [1] we also get a better 
empirical upper bound on the true ground state energy 8,g{N) than £.g{N) from 
the remaining data of the computer-generated list {£^(A^); A^ = 2, 3, A^*}: 

Note that (HI) is always true, but it only leads to a better empirical upper bound 
on £g(A^) than S,g{N) when ^g{N) fails a test. Note also that the computer- 
experimental data supply empirical upper bounds to the actual ground state 
energies within the numerical accuracyp i.e. we always have £^(A^) > E,g{N). 

Subjecting some publically available lists of computer-generated data of 
putative ground state energies £.g{N) for various A^-body systems to the above 
tests has yielded several dozen instances of £.g{N) which failed one of these 
tests; i.e., for those A^ we conclude that £g(A^) > ^g{N) strictly, and we get an 
improved empirical upper bound on ^g{N) through (jlj). The surveyed A^-body 
systems are A^ point charges with D-dimensional Coulomb pair interactions 
(D = 2; 3), placed either on the unit 2-sphere or on a 2-torus (variably known 
as (elliptic) Fekete, Thomson, and Riesz problems). The analysis of the sphere 
data is reported in the next section; for the torus data see section 3. 

^It is assumed here that the algorithms have computed the energy of some configuration. 
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2 Many point charges on the 2-sphere 

Finding the minimum energy configuration(s) of point charges placed on 
the unit 2 sphere is a beautiful, intriguingly rich, and hard mathemati- 
cal problem which in addition is relevant to many fields of science; see the 
survey articles |ErHo97j . |SaKu97] . |AtSu03j . and |HaSa04] . and the website 
|Wom09] ■ One can either interpret §^ as a two-dimensional "physical space" 
in its own right with f/§2(qr^, qg) given by the D = 2-dimensional Coulomb pair 
interaction —ln\q^ — q2\, where \qi — q2\ is the cordal distance on inciden- 
tally, the cordal distance on coincides with the three-dimensional Euclidean 
distance between two points and ^2 C M^, but the embedding can be 

avoided in the discussion. Or, one can interpret §^ C MP for D > 2 as a proper 
submanifold, with Us2{qi, ^2) given by the D-dimensional Coulomb pair inter- 
action |q^ — 92^ '^! where |q^ — 92I is the D-dimensional Euclidean distance 
between two points and ^2 C MP. For small N the ground state 

configuration (a.k.a. an A^-tuple of Fekete points) can easily be characterized 
explicitly!! and the asymptotic large A^-dependence of ^g{N) can be, and to 
some extent has been |RSZ94j determined analytically without seeking the 
exact Fekete pointsl^l but in general the problem defies analytical treatment. 
Computer experiments (e.g. |RSZ94j . |RSZ95j . |Aetal97j . |Petal97j . |ErHo97j . 
|BCNT02j . |BCNT06j . |Betal07j ) help finding candidates for the minimizing 
configuration and in any event yield empirical upper bounds £^(iV) on the 
ground state energy 8,g{N). But even computers are soon overwhelmed be- 
cause the number of local minimum energy configurations which are not global 
seems to grow exponentially with N |ErHo97j so that a computer algorithm is 
more and more likely to find one of these non-global minima. Indeed, our tests 
have successfully detected a couple dozen non-global minimum energy values 
because their pair-specific value surpassed an ensuing empirical pair-specific 
energy value in some computer-generated list. 

2.1 Two-dimensional Coulomb interactions 

Tables of computer-experimental ground state energies {8,g{N)] N = 2,3, 4, ...} 
for point charges on with two-dimensional Coulomb pair interaction 
Us2{qi, ^2) = — In IQi — Q2I can be found on the website |BCM] . more precisely: 

|http: / / physics.syr.edu/condensedmatter/thomson/ shells/ ... 

...shelltable.php?topology=sphere&potential=0&start=0&end=5000, 
and in |RSZ95j . which contains many references to earlier studies. 

^For D = 2&3 the first dozen minimizers are discussed in the delightful article |AtSu03| . 

''it "suffices" to know that for large N the Voronoi cells around the charges are mostly 
hexagons of a certain size; see |SaKu97j for an enlightening discussion. We also remark that 
the leading order term in the asymptotics can be determined with the help of a very general 
variational principle, see [KiSp99| . 
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A value of ^^{N) is listed in |RSZ95] for 199 consecutive values of N 
starting at = 2; one also finds there three additional ^g{N) values for 
N which correspond to the icosahedral group. On p. 116 of their article, the 
authors describe various careful tests which they have undertaken to increase 
the chances of their £^(A^) actually being 8.g{N). In addition to their tests I 
now have subjected their list of data {8.g{N)]N = 2,3, ...,200} to the n = 1 
monotonicity test (which, apparently, is not among the tests which the authors 
of |RSZ95] have used). Their data passed the n = 1 test (there is no point, 
then, to test with n > 1). 

While all tables in |RSZ95j are of course permanent, those at |BCMj . ac- 
cording to these authors, are interactive and are updated whenever some user 
finds a new and better pair (A^, 8,g{N)). Therefore it is mandatory to also give 
the information on which day one downloaded the data from [BCM] for study, 
which for the computer-experimental data on the D = 2 Coulomb ground 
state energies I did on Feb. 21, 2009; yet, on Feb. 24, 2009 the crucial data 
I will be talking about below were still the same. At the time of writing of 
this paper, this data hst {£^(A^); A^ G {3, 5000}} has (had) plenty of gaps, 
i.e. putative ground state energies 8,g{N) are absent for many values of A^. 
While the absence of the case A^ = 2 is not a problem, as £^(2) is exactly 
known, the first real gap is the absence of any £^(11), any £^(19), and soon 
the gaps become larger and larger. The larger the gaps, the less likely one is 
to detect non-globally minimizing &g{N) with the monotonicity tests, yet two 
data points "got caught in this net": £^(97) = -891.653265231 > £g(97), for 
e^(100)-e^(97) = -0.013678811 < 0; and £^(2000) = -386, 187.080630499 > 
£,,(2000), for £^(4212) - £^(2000) = -0.000503199 < 0. Notice that while 
there is only a gap of 2 missing data between £^(97) and £^(100), the gap 
between £^(2000) and £^(4212) is a whopping 2211. Since the non-globally 
minimizing £^(2000) was detected by the n = 2212 test, it must be quite 
far away from the actual ground state energy £^(2000). By (jl]), and using 
£^•(4212) = -1, 722, 205.927290610 from [BUM] , we find the upper bound, 

£9(2000) < miff £^(4212) = -388, 198.8687. (5) 

Similarly, with £^(100) = -1,083.376338235 from jBCM] . withn = 3yields 
£3(97) < ^^£^(100) = -1,019.030349; however, this upper bound is cer- 
tainly beaten by £^(97) = -1, 022.023977757 in |RSZ95j . Incidentally, several 
other non-globally minimizing data £^(A^) in the list at |BCM] which actually 
pass the monotonicity tests can be detected by simply comparing with the list 
in |RSZ95j (at the time of writing, the non-global data at |BCM] are £^(12), 
£^(36), £g(60), £^(87), £^(96), £^(100)); yet, no such comparison is possible 
for £^(2000) as the data reported in |RSZ95j do not go beyond A^ = 282. 

We end this subsection with an illustration of the monotonicity of the 
sequence A^ 1— > Bg{N) by plotting the monotonically increasing sequence A^ 1-^ 
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£g{N) > eg{N) using the data in |RSZ95j . with "=" for < 15 |AtSu03j . 




+ 

-0.3- 



■ + 
-0,4- 

Computer-experimental pair-specific ground state energies Eg (N) for N point charges 
on with logr~^ Coulomb interactions, using the data of RSZ95 (crosses). Shown 
also in this plot is the large N asymptotic two-term approximation to eg{N) (dots). 



Aheady for > 50 this empirical curve seems to agree to within less than 1% 
absolute error with the large asymptotics of the actual sequence N £g{N), 
given by eg{N) x ^ In | — |iV~^ In iV + 0{N~^) and obtained through dividing 
by N{N-1) the large-iV asymptotic formula for N ^ 8,g{N) |RSZ94[[RSZ95] . 

Eg{N) = aN^ + bN\nN + 0{N) (6) 

with a = I In I and b = —\. In |RSZ94j it is also conjectured that, actually, 

Eg{N) = aN^ + bN\nN + cN + d\nN + 0(1), (7) 

and rigorous upper and lower bounds on c are given there. Smale's 7*^ problem 
for the 21** century asks for an algorithm which upon input returns an Ap- 
point configuration on for which 8,g{N) does not deviate from ^g{N) by 
more than the fourth term in the conjectured expansion ([7]) (possibly up to a 
different coefficient d'), and which does so in polynomial time; see |Sma98j . 



2.2 Three-dimensional Coulomb interactions 

We next discuss the data for A^ point charges on the sphere with three- 
dimensional Coulomb interactions f/§2(q^,qr2) = Vl^i ~ since |Why52] 
referred to as "Thomson's problem" even though Thomson's "plum pudding 
model of the atom" [Tho04 j is not quite the same problem. Be that as it 
may, there seem to be many more studies of this Thomson problem than of 
its variants with other pair interactions and geometries. I perused a sample of 
those studies and eventually found data which failed a monotonicity test. 

Starting at A^ = 2, in |ErHo97] one finds 111, and in |RSZ95] 199 con- 
secutively computed putative ground state energies £^(A^) for the Thomson 
problem (as just defined). Since the authors in |RSZ95j mention that within 
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numerical precision their data agree with those on Sloan's home page |HSS94] 
(actually, at the time: its predecessor via netlib.att.com), I opted for analyzing 
the data on |HSS94] which are stored (much) more userfriendly than those in 
|RSZ95j and |ErHo97] . even though for only 129 consecutive values of some 
value £g(iV) is listed, starting with = 4, plus a handful of other values for 
A^s up to = 282. All data at |HSS94j passed the n = 1 monotonicity test. 

A wonderful treasure of data for the Thomson problem is listed on the 
website jBCM] , more precisely: 

|http: / / physics.syr.edu / condensedmatter / thomson / shells / '. . . 

...shelltable.php?topology=sphere&potential=l&start=0&end=5000. 
I downloaded the data on Feb. 09, 2009, and on Feb. 24 the relevant data were 
still present. Starting with A^ = 3, one finds computer-experimental ground 
state energies £^(A^) for an amazing 1488 consecutive values of A^, the first 
gap occurring at A^ = 1491. In addition one finds almost 300 wider spaced 
^g{N) up to A^ = 5000. Quite remarkably, all the consecutive values up to 
A^ = 1490 passed the monotonicity test, and so did an additional 70 values of 
^g{N) which are listed with gaps up to A^ = 1800. The computer-experimental 
£^■(1801) = 1,579,605.0292504800 is the first in the list which failed a mono- 
tonicity test, in fact the n = 1 test, with 5^(1802) -e^(1801) = -0.0000044325. 
Of course, there is also a better upper bound on £^(1801) via (jl]), but it 
is not much of an improvement. The story is quite different for the next 
failling, which is £^(2002) = 2,004,888.5938241700, which failed the lO*'' and 
the 20*^* monotonicity tests, with £^(2012) - e^(2002) = -0.0125431412 and 
£^(2022) - £^(2002) = -0.012518560. This time we obtain two upper bounds 
from (jlj) (for n = 10 and n = 20), and the better one is for n = 10, viz. 

£,(2002) < ffffi £^(2012) = 1, 954, 640.745 (8) 

which is also a considerable improvement over the listed £^(2002). Not sur- 
prisingly, in this long list of computer-experimental data our harvest is richer 
than in the previous log list. An additional nine valus of &g{N) failed one or 
more of the monotonicity tests (indicated by (n = x) behind £^(A^)), namely 

£^(2531) = 3, 204, 550.3074368200 (n = 1; 19) 

£^(2561) = 3, 207, 772.6856023400 {n = 1) 

£^(3362) = 5, 543, 845.2403717400 {n = 2) 

£^(3480) = 5, 941, 792.7610333500 (n = 1; 2) 

£^(3663) = 6, 586, 476.2826423300 {n = 9) 

£^(3720) = 6, 793, 857.9289983900 (n = 1) 

£^(4000) = 7, 860, 293.8236758000 {n = 2) 

£^(4260) = 8, 920, 193.5261720900 {n = 2) 

£^(4620) = 10, 498, 739.0438109000 {n = 2; 4). 
We leave it to the interested reader to use (jl]) to compute upper bounds on 
^g{N) for the A^ showing in this table and the pertinent £^(A^-|-n) at [BCM] . 
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All the data in |Aetal97l[Petal97] pass the monotonicity test, but these data 
are few and far between so that it is actually quite unlikely for any of them to 
fail a monotonicity test. In fact, [ Aetal97j in their table list only five different 
N e {2,...,2500}(!), with two different 8-giN) each; |Petal97] point out that 
they found a lower energy £^(2472) than did [Aetal97j . Also [BCMj provided 
a putative value for = 2472, and their £^(2472) = 2, 987, 485.953(...) is 
even lower than the one in |Petal97j . which is £^(2472) = 2,987,486.132. 
Incidentally, there is also a moral here, to be told in the last section. 

To illustrate the monotonicity of the sequence i-^ ^g(^) I plot the mono- 
tonically increasing sequence A^ i— >• Sg{N) > eg{N), this time using the data 
in |HSS94j . For A^ < 15 also these two sequences are identical [AtSuOS] . 




Computer-experimental pair-specific ground state energies eg{N) for point charges on 
with Coulomb interactions, using the data of [HSS94 (crosses). Also shown are the 
leading three terms in the conjectured asymptotic large A' approximation to eg{N) (dots). 

Also shown in this diagram is the partially conjectured large A^ asymptotics 
eg{N) X 1 - 0.55305(...)A^-i/2 ^ i^-i + 0(iV-3/2), obtained through dividing 

by A^(A^ — 1) the following large- A^ asymptotic formula for A^ \—>- 8,g{N), 



^g{N) = aN^ + bN^/^ + cN + dN^/^ + e + 0{N-^ 



(9) 



where a = 1/2 is the only rigorously proven coefficient |RSZ94] . |KuSa98j . 
|KiSp99| , while it is conjectured |RSZ94j that c = = e and 
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whereas d is estimated numerically in |RSZ95] . 



3 Many point charges on a 2-torus 



(10) 



For the perhaps most prominent non-spherical topology, the 2-torus C I 
we found data lists at jBCMj . However, these lists are clearly preliminary. 
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3.1 Two-dimensional Coulomb interactions 



Curiously, at |BCM] putative ground state energies ^g{N) are only listed for 
four different values of A^, and all of them are obtained with different aspect 
ratios of the tori. Hence, no meaningful monotonicity test can be applied. 

A visualization of a putative ground state configuration of point charges 
with logarithmic Coulomb interactions on a 2-torus can be found in |HaSa04] 
and on the cover of that issue of the Notices, and also at |Wom09j . 



3.2 Three-dimensional Coulomb interactions 

At |BCM] one finds about 50 data of putative ground state energies E>g{N) 
for the aspect ratio 1.414, which are computed for sparcely placed A^ up to 
A^ = 5000. For sparsely placed data one would expect it to be less likely to 
find some which fail a monotonicity test. However, the data £g(A^) for nine A^ 
failed a monotonicity test for one or more n, namely: 

£^(15) = 81.390479174 (n = 5; 6; 9; 11; 12) 

£^(30) = 331.088832684 (n = 2; 5; 6; 7; 10) 

£^(113) = 5, 370.892624565(^ = 4; 7; 38; 49; 78) 

£^(262) = 28, 287.128667479 (n = 38; 238; 290) 

6^(360) = 53, 857.158956562 {n = 140; 192; 200; 694) 

£^(396) = 66, 660.796433247 {n = 104; 156; 164; 598) 

£^(1000) = 418, 396.928796506 (n = 363; 4000) 

£^(1363) = 707, 154.008010865 [n = 3637) 

£^(3500) = 5, 174,438.587013800 (n = 1500). 
Again, we leave it to the interested reader to use (jl]) and the above table to 
compute upper bounds on S,g{N) from the pertinent data list at |BCM] . 



4 Some variations on the theme 

Two- and three-dimensional Coulomb pair interactions and 2-sphere and 2- 
torus domains are merely the most prominent examples of pair interactions 
U\{q^, Qj) and d-dimensional domains A to which Proposition[T]applies. Other, 
though physically less important, examples of interactions are D-dimensional 
Coulomb interactions with D > 3, and more generally the so-called Riesz in- 
teractions, computed with uj^\q^, Qj) = —sign {s)\q^ — qj\'^, for any real s < 2; 
see, e.g. |KuSa98j . |HaSa04j . |HaSa05j . |KSS07j . |Betal07j . The logarithmic 
Coulomb interaction is usually considered to be the special case s = 0, in the 
sense that lim^j^o ^""^(Iqj — qj\~'^ — 1) = — In Iq, — q^l; For s = 1 the Riesz en- 
ergy gives U\{qi, q^) = —\qi — qj\, in which case Proposition [T] may shed some 
new light on the question of the maximum average pair-wise distance of points 
in A and related problems; beside the cited general survey articles, also see 
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|Bec84j . Proposition [T] applies also to other bounded domains, in particular 
curves! 1 should emphasize that the logarithmic interactions between charges 
constrained to (planar) curves can be studied in quite some detail with complex 
variable techniques, see |Ketal04j . Proposition [1] can easily be generalized to 
unbounded domains, with lower semi-continuity replaced by another appropri- 
ate condition guaranteeing minimizing configurations for all A^. A physically 
important example is A = with U^i^q^.q^) = \qi ~ qjl^^"^ — \qi — Qjl^^y 
which has minimizing A^-body configurations for each A^, known as Lennard- 
Jones clusters, see |AtSu03] for a recent survey. At the expense of replacing 
the minimum by an infimum, pair interactions which are merely bounded be- 
low can be handled also, but minimizing sequences which don't converge to a 
minimizing configuration are perhaps less interesting. 

5 Summary 

The main purpose of this article is to draw attention to the monotonicity 
tests implied by Proposition [1] and to emphasize the ups and downs of this 
monotonicity test family, not to report on an exhaustive series of such tests 
covering all available data. In fact, having demonstrated the utility of these 
tests, it is much more efficient when they are being directly implemented in 
the computer experiments rather than being run by a third person afterwards. 

1 still owe the reader the moral announced earlier. Originally 1 had analyzed 
the data of computer-experimental ground state energies ^g{N) reported in 
[Aetal97j which, divided by A^(A^ — 1), arranged themselves monotonically 
increasing when plotted vs. A^; in fact, this prompted me to conjecture and 
then prove Proposition[Tl But, as we saw, those data are not the correct ground 
state energies. With hindsight 1 was quite lucky, for such widely spaced data 
almost inevitably form a pair-specific monotonic sequence. Had 1 hit upon a 
more closely spaced list of non-optimal data which would not have been pair- 
specifically monotonic, 1 may not have conjectured the monotonicity in the 
first place! 

Acknowledgment. This paper was written with support from the NSF under 
grant DMS-0807705. Any opinions expressed in this paper are entirely those 
of the author and not those of the NSF. 1 like to thank my colleague A. Shadi 
Tahvildar-Zadeh for going through the painstaking process of verifying that 
those data which failed a monotonicity test were really there and not made up 
by meQ 1 also thank Jozsef Beck and Doron Zeilberger for their comments on 
Proposition 1, and Xiaoping Peng for helping me with the data processing. 

^ Just in case you take a look and don't find those data anymore which failed a monotonic- 
ity test — there were some that did! Incidentally, I myself plan to try to improve on some 
data which failed a monotonicity test after registering on the interactive website |BCMj . 
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